COMBINATORIC RESULTS FOR GRAPHICAL ENUMERATION AND THE HIGHER 

CATALAN NUMBERS 



VIRGIL U. PIERCE 

Abstract. We summarize some combinatoric problems solved by the higher Catalan numbers. These 
problems are generalizations of the combinatoric problems solved by the Catalan numbers. The generating 
function of the higher Catalan numbers appeared recently as an auxiliary function in enumerating maps 
and explicit computations of the asymptotic expansion of the partition function of random matrices in the 
unitary ensemble case. We give combinatoric proofs of the formulas for the number of genus and genus 1 
maps. 



1. Higher Catalan Numbers 

The Catalan numbers solve a number of classical combinatoric problems such as the "Euler Polygon 
Division Problem" : how many ways are there to divide a marked polygon with j + 2 sides into triangles 
using edges and diagonals [31 [71 [HI [121 [IS] (see figure [T|) . 




Figure 1 . A polygon with 4 + 2 = 6 sides divided into 4 triangles using edges and diagonals 



They count the number of right-left paths along a 1-Dimensional integer lattice which stay to the right of 
and go from to in 2j steps; equivalently they count Dyck paths from (0,0) to (2j,0) [T1 l4l ITSl 118]. 
They count the number of ways for 2j customers to line up, with j customers having only a 2-dollar bill and 
j customers having only a 1-dollar bill, to purchase l-doUar widgets, so that each customer receives exact 
change. They count the number of non-crossing handshakes possible across a round table between n people 

m- 

In this paper we will explore a generalization of the Catalan numbers, the higher Catalan numbers. 
We will show that this generalization solves enumerative problems that are natural generalizations of the 
problems solved by the Catalan numbers. We will then highlight their appearance in recent results on map 
enumeration problems. 



Let 



(1.1) z(,s) = l + 5^C 

be defined implicitly as the solution to the algebraic equation 

(1.2) szisY - z{s) + 1 ^ 
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which is regular at s = 0. This generating function was presented in [171 [20]. It is known that the (j'''' are 
(see [Tl[20]) 

(1.3) cr' 

The numbers p.3p are called the higher Catalan numbers and solve a number of combinatoric problems. 

The C}^ are the Catalan numbers. The higher Catalan numbers count the number of ways to divide a 
marked {{i' — l)j + 2)-sided polygon into (i^ + l)-sided polygons using edges and diagonals as in figure [5] 



Figure 2. A polygon with 2 • 5 + 2 = 12 sides divided in to 5 squares, i.e. = 3, using 
edges and diagonals 



They count the number of paths along the 1-Dimensional integer lattice which move to the right 1 unit 
per step, to the left (j/ — 1) units per step, stay to the right of and go from to in vj steps [TS]. They 
count the number of higher Dyck paths: paths which go from (0, 0) to {lyj, 0) along the 2-D integer lattice 
with steps (1, 1) and (1, — (j/ — 1)) which stay above the x-axis. They count the number of ways for (ly — l)j 
customers with l-doUar bills and j customers with i^-doUar bills to line up to buy l-doUar widgets so that 
each customer receives correct change. They count the number of non-crossing handshakes possible across a 
round table between n i/-handed beings. 

Remark 1.1. More generally we call generalized Dyck paths walks on the 2-D integer lattice restricted to the 
upper half plane with steps coming from a set of integer vectors. The generating function of the number of 
generalized Dyck paths satisfies a system of algebraic equations [13] . 

In addition to these combinatoric problems the coefficients of z(s) satisfy the recursion relation 

(1.4) cl^)= c^'c^'---c^- 

jl+j-l^ hj,^=j-l 

These numbers appeared in counting labeled maps using the partition function of random matrices [IT]. 
We found, in solving the map enumeration problem for some simple cases, that the generating function could 
be written as a function of the generating function for the higher Catalan numbers. 

A map D on a compact, oriented connected surface AT is a pair D = {K{D), [t]) where 

(1) K{D) is a connected 1-complex; 

(2) [l] is an isotopical class of inclusions l : K{D) — > X; 

(3) the complement of K{D) in a; is a disjoint union of open cells (faces); 

(4) the complement of Kq{D) (vertices) in K{D) is a disjoint union of open segments (edges). 

Maps can be thought of as ribbon graphs (or fattened graphs) embedded on an oriented surface [H [THl [H] . 

Define eg{t) to be the generating function of the number of connected maps of degree 2v embedded in a 
genus g surface, in the following sense 

e.w=f:4^^H.)^f, 

where ti'g^^j) is the number of planar maps with j vertices of degree 2v embedded in a genus g surface. 
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Define the class of functions: iterated integrals of rational functions of z, or iir of 2, to be functions which 
are found by taking a finite number of anti-derivatives of a rational function of z. The particular subclass 
we are going to be concerned with are those izr's which have singularities only at z = 0, 1, or y / {v — 1). One 
may check that this subclass of functions is closed under anti-differentiation. 

To illustrate the results of [TT] we present the explicit result we found for 5 = 0. Let 



Cy — 2v 



2v-l 



Theorem 1.2 (Ercolani- McLaughlin-Pierce) 



eo(-t) = Hi,{z{cJ) - l){z{ci,t) -rt,) + - \og{z{cJ)) , 



where 



and 



4i/(i/-h 1) 
3(i^ + l) 

= — • 

V — 1 

More generally we find that eg{—t) is an iir function of z{cyt) with singularities only possible at z — Q, 1, or 
vl{v-\). 

We will now outline the approach we used to prove this result. Define the partition function of random 
matrices 



(1.5) 



ZN{ti,t) = J cxp ^-iVTr (^tiM + ^M^ + tM^''^ ^ dM, 



where the integral is taken over the space oi N x N Hcrmitian matrices, and dM is the product of Lebesgue 
measures on the independent variables in M: 



dM = 



[]d(Rc(Afy))rf(Im(My)) 



Let 



/bi ai 



L 



ai 62 

02 



02 
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Y[dM„ 

. i 

J 



The Toda lattice hierarchy is given by the differential equations 

dL 
dt-i 



where Bj = {L^)+ — (L^)_, and where {■)± indicates projection onto the upper (resp. lower) triangular parts. 
This hierarchy is an integrable hierarchy, possessing a complete family of independent commuting integrals. 
Solutions are generated as logarithmic derivatives of so called r-functions. The partition functions (|1.5p are 
these r-functions. 

Using this fact we showed that the leading order (as N 00) of a second logarithmic derivative of the 
partition function was the generating function of the higher Catalan numbers. 



Theorem 1.3 (Ercolani- McLaughlin-Pierce). 
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zi-ca) + o 



Ar2 



27V2 dt\ 



log 



^jv(0,0) 



ZN{N^''^ti,N^-''t) 



ti=0 



The paper [TT] presents two independent proofs of Theorem 11.21 in the case when 5 = 0: one is a direct 
calculation of the generalization of the method of steepest descent for computing the leading order of asymp- 
totics of matrix integrals, the other uses the fact that the partition function (jl.Sp is a r-function of the Toda 
lattice hierarchy. Theorem 11.31 is found in either proof of Theorem II .21 in the g = case. The general result 
from this paper is that eg{t) is explicitly computable as a function of z{s). Therefore the higher Catalan 
numbers play a central role in the combinatorics of maps. 

1.1. Results. We will summarize our results for enumerations of genus and genus 1 maps here. The 
generating function of the number of genus maps is 



— -tj 

while the generating function for genus 1 maps is 



J- 



We will give a combinatoric proof of the theorem from fllj: 
Theorem 1.4 (Ercolani-McLaughlin-Pierce) . The Taylor coefficients of bq are 

S2u) _ ciiyj - 1)! 



and those of ei are 

In section [5] we will review the combinatoric arguments showing that the higher Catalan numbers are 
generated by z(s). Both the path counting problems analyzed in section [2] and the queuing problem men- 
tioned are from Yaglom and Yaglom |2QJ, and Sato and Cong ,18J. The path counting problems have been 
cosmetically changed, but the method of solution is still the same. In section[3]we give the origin of theorems 
like 11.21 and our general result. In section |4] we give a combinatoric calculation starting from Theorem 11.21 
which culminates in a formula for the number of planar maps with j vertices of degree 2v. In section [S] we 
carry out a similar calculation for counting genus 1 maps, giving a proof of Theorcm ll.4l 

2. Calculation of the Higher Catalan Numbers 

The higher Catalan numbers satisfy a recursion relation which will tie together four types of combinatoric 
problems: the Taylor coefficient problem (jl.ip . the polygon division problem, and two problems we are 
about to introduce: a path counting problem and a queuing problem. The path counting problem is directly 
solvable by enumeration. This is then connected to the queuing problem and we show that both satisfy the 
same recursion (|1.4p as c}^^ and the polygon division problem. 

We give here a modification of the calculation outlined in Yaglom and Yaglom [50], and Sato and Cong 
[T5] . In this section we will show that the higher Catalan numbers, 

count the number of right 1, left (j/ — 1) walks on the 1-Dimensional lattice which go from to in lyj steps 
which stay to the right of 0. 

Note that this 1-Dimensional counting problem is identical to counting the number of (1, 1) and (1, — (j/ — 
1)) paths on the 2-Dimensional Integer Lattice from A = (0,0) to B = (i^j, 0) which stay above y = 0, we 
call these higher Dyck paths. 

As an example of how to count such paths we have the lemma 
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Lemma 2.1. The number of (1,1) and — 1)) paths from P = (pi,P2) to Q — (<7i,<72) ( assume 

Pi < Qi) on the 2-Dimensional integer lattice is 

AT f^^^ Pi 

where n, u, and d solve the system of equations 

Pi + u + d = qi , 
P2 + u ~ {ly - l)d = q2. 
If no integer solutions exist then there is no such path connecting these points. 

To count the number of(l,l), (1,— (^— 1)) paths from A to B which stay above the x-axis we will count 
the number ol such paths which pass below the x-axis. A path of this type going below the x-axis will go 
through one of the points 

Dk = (i^fc- 1,-1). 

Let A' = {0,-iy). 

Lemma 2.2. The number of (1, 1), (1, — (i/ — 1)) paths from A to B which pass below the x-axis is 



{v-i)NA'B = {y-i) 

Begin the proof of this lemma by computing 

Nad, - 

and 



vk-l 
k 



Now consider the number of paths from A to B which go through Di 

Nad,Nd,b = {y- 1)Na'd,Nd,b, 

which is (z^ — 1) times the number of paths from A' to B through Di. Then inductively assume that the 
number of paths from A to B through Di is (z^ — 1) times the number of paths from A' to B through Di for 
I < k. Consider the number of paths from A to _B which go through Dk but not Di for I < k, 

Nad,Nd,b - Nad,^,Nd,_,d,Nd,b - ■ ■ ■ = {,y - 1)Na'd,Nd,b ' - 1)Na'd,^,Nd,_,d,Nd,b 

□ 

The number of (1,1), 1)) paths from A to B which stay above the x-axis is 

(7)-<''-"G-^)4G-^, 

which is the j'th j/-higher Catalan number. 

Consider the queuing problem: How many ways are there for (j/ — l)j customers with 1-dollar bills and 
j customers with j/-dollar bills to line up to buy l-doUar widgets so that each customer receives correct 
change. In the path counting problem, in order for a step to the left to be placed on the path and stay to 
the right of there must be (i^ — 1) corresponding right steps that have come before on the path. In the 
queuing problem, in order for a customer with a z/-dollar bill to receive correct change there must be (i/ — 1) 
corresponding customers with l-doUar bills in line ahead of the i/-dollar bill. 

To see that the numbers (|2.ip agree with the coefficients (|1.3p of z{t) we argue that the number of right 
1, left {ly— 1) paths which stay to the right of satisfy the recursion relation (|1.4p . As argued above, to each 
left step in the path there must correspond (j/ — 1) steps to the right. The path begins with a step to the 
right. This step corresponds, together with {v — 2) other right steps, to a left step later on. This collection 
of p steps divides the entire vj steps into sub-paths of the type: steps right by 1, left by {v — 1), from k 
to k which stay to the right of k. Therefore the recursion relation p.4p is satisfied. 

Now we argue that the number of ways to divide a marked ((j^ — l)j + 2)-sided polygon into (y ^ l)-sided 
polygons using edges and diagonals satisfies the recursion relation (|1.4p . Take the sub-polygon which has the 
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marked outside edge as an edge. This sub-polygon divides the {{ly — l)j + 2)-sided polygon into v polygons, 
therefore the recursion relation (|1.4[) is satisfied by this problem as well. Therefore these numbers are also 
the higher Catalan numbers. 

3. The Enumeration of Planar Maps 

Our interest in these combinatoric problems began with the function z{s) defined by (|l.ip . This function 
appeared as an auxiliary function in explicit calculations of the asymptotic expansion of the partition function 
of large random matrices. It is interesting that a generating function as classical as z{s) appeared in this 
setting as the terms of this asymptotic expansion are themselves generating functions for a combinatoric 
problem. The Taylor coefficients of the terms in the asymptotic expansion enumerate maps on genus g 
surfaces. 

In this section we highlight the information contained in the asymptotic expansion of the partition function 
(jl.Sp for large matrices. We have detailed data about the higher genus problems [TT], however here we will 
concentrate on the genus zero (or planar) setting. 

Recall that the partition function of random matrices we consider is (jl.5[) : 

ZN{t) = J cxp (^-NTi + tM^"^^ dM, 

where now we take ti = 0. Ercolani and McLaughlin lOj showed that \og{ZN{t)/ZN{0)) possesses an 
asymptotic expansion as ^ oo inside a non-trivial t domain, 

/ ^ 1 , ( ZN{t)\ , , 1 , , 1 , , 

(3.1) ^ log j ^ eo(t) + -e.it) + -e.(0 + . . . 

where eg{t) is analytic in a neighborhood of t = and is a counting function for genus g maps with 2v- 
degree vertices The idea of counting maps in this way originates with random matrix models of quantum 
field theories [21I61II9]. 

In this paper we will concentrate our attention on the genus and genus 1 terms in p.ip . The genus 
term is 

where the h^q^^ (j) are the number of genus maps with j vertices of degree 21/. These numbers are calculated 
in |llj by explicit computation of the contour integral 

(3.2) ^(^■^)^j)^jL<j;t-=-ieo{-t)dt, 

where the contour encircles t — and cq is given in Theorem 11.21 To evaluate (|3.2p one rewrites it as an 
integral with respect to z = z{ci,t). The result is the first part of Theorem 1 1.41 

4. Computation of the Taylor Coefficients of eo 

We will now compute the Taylor coefficients of eo from combinatoric arguments only. In section [2] we used 
combinatoric techniques to compute the Taylor coefficients of z{s). Theorem 11.21 gives eo(— i) as a function 
of z{s). The calculation of the Taylor coefficients of eo(— t) can be done from this theorem utilizing contour 
integration and a clever change of variables. Our goal in this note is to complete this calculation using 
only combinatoric arguments. To that end, we need to derive formulas for the coefficients of log(z(s)) and 
{z{s)-l)\ 

4.1. Taylor CoefRcients of powers of {z{t) — 1) and log(z(i)). The coefficients of (z(s) — 1)* are given 

by 

The flf'^^ count the number of right 1, left {v — 1) paths from to in vj steps which stay to the right of 
and return to at least i — 1 times in between the ends of the path. This is equivalent to the queuing 
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problem: how many ways can j customers with I'-dollar bills and — l)j customers with 1-dollar bills form 
i lines to buy 1 dollar widgets so that each customer receives exact change. 

This multi-line queuing problem is equivalent to: how many ways can j — i customers with i/-dollar bills 
and (i^ — l)j + 1 customers with 1-dollar bills form one line to buy 1 dollar widgets so that each customer 
receives exact change. This is done by shifting the lines together. 

First, we notice that at the back of the first line is a z^-doUar bill. We shift the problem in the following 
way: we change the i^-doUar bill at the back of first line into a 1-dollar bill, then adjoin the second line to 
the back of the first. The first line now has j — 1 customers with ^-dollar bills and {ly — l)j + 1 customers 
with 1-dollar bills (a i^-doUar bill has been changed into a 1-dollar bill) . We are now in the case of i — 1 lines 
and may shift again. After i — 1 repetitions of this process we have a single line arranged in such a way that 
the last customer has a z^-doUar bill, we remove this customer from the process. 

To show that the two counting problems are equivalent we show that the one-line queuing problem above 
can be transformed back into the multi-line problem. We begin with a line of j — « customers with j/-dollar 
bills and {f — l)j + 1 customers with 1-dollar bills in a single line to buy 1 dollar widgets so that each 
customer will receive exact change. First add a j/-dollar bill to the back of the line. Starting from the back 
of the line and moving forward we will form sets of blocks of v customers. Each block has a customer with 
a j/-dollar bill and i/ — 1 customers with 1-dollar bills. When we reach a customer with a 1-dollar bill for 
which there is no corresponding i^-doUar bill behind, we have completed a line, and we remove it from the 
first line. We then turn the 1-dollar bill at the back of the first line into a i^-doUar bill and proceed as before. 
Thus transforming us to the multi-line queuing problem. 

The path counting problem which corresponds to this shifted queuing problem is: how many right 1, left 
(i^ — 1) paths, from to — l)i-|-l, in i/j — z-|-l steps, are there which stay to the right of 0. This is equivalent 
to the problem: count the number of (1, 1), (1, —{v — 1)) paths from (0, 0) to {vj — i + 1, {ly — l)i + 1), which 
stay above zero. These paths are counted in the same way as the ones in section 2: 

(4-1) % ' ]> I- 

The coefficients of log(z(s)) are computed directly using the Taylor expansion of the logarithm function: 

iogMs))^E^r^^' = E(-ir ^'°^'^" ^ - 

We find that the coefficients are 



3 \3 - V i V J - 1 



4.2. Enumeration of planar maps. Combine the results (|4.ip and (|4.2p with Theorem 1 1.2 

eo(-i) = fi{z{cj) - l){z{ci,t) - r) + i \og{z{c,yt)), 



and we find that 

''Hj) = 



"((^-1)^ + 2)!- 
This proves the first part of Theorem 11.41 

4.3. Additional enumerative results for functions of z{s). Formulas similar and equivalent to (|4.ip 
have been derived for other functions of z(s) ^M,- 

Theorem 4.1 (Polya-Szego-Gould ). Let z{s) be the generating function for the Higher Catalan Numbers. 
Then 

oo 

z"=^A,(a,^.)s^ 
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with 

(4.3) A^{a,v) ^ -———[ 

In addition 



l)z ^ 



I' — (ly — l)z ^-^ V 7 



One can check that formula (|4.3p is equivalent to (|4.ip . 

5. The enumeration of genus one maps 

In the paper [TT] the authors show that there is a general construction of the functions eg(tiv) in the 
asymptotic expansion p.ip as functions of z{cii,tii,') . For example, we find that in the case of genus one 
maps 

Theorem 5.1 (Ercolani-McLaughlin-Pierce ). In the neighborhood of analyticity, the second term in the 
asymptotic expansion Ii3.1\) is given by 

(5.1) ei(-t) = log l)2(c2.t)) , 
where we choose the principal branch of the logarithm. 

Our goal is to derive a formula for the Taylor coefficients of ei (— ^2;^) for general ly. This is a straightforward 
calculation: expand formula (|5.1|1 as a powers series in [z{c2vt) ~ 1); 

(5.2) ei(-t) ^ J^^^^ (^(^2-^*) - • 

k=l 

Then insert the power series of (z — 1)''' into (|5.2p : 

fc=l k=l j=k J \J / 

where the rj^'^''^'^ were computed in (|4.ip . Switching the order of summation in (|5.3p we find that 



(5.4) {j) = ^{1^-1, , 



This proves the second part of Theorem 11.41 

One is tempted to try to use the last part of Theorem 14.11 to find a more concise formula for n'^'^^j). 
However using the relation 



u — {v — l)z ' 



derived from the implicit definition of z by formula (|1.2p . we see that the derivative of ei(— t) with respect 
to t is 

(5.5) eU-i)-^^-^) 



12 l)z{ca)f 



There does not seem to be a more concise formula for the coefficients of (j5.5p than that found directly from 
j times K^i^^ in (|5.4p . 



6. Conclusion 



The higher Catalan numbers appeared naturally in the study of the enumeration of maps embedded 
on Riemann surfaces. We have gathered here a number of interesting facts about these numbers and the 
combinatoric problems related to them. It would be interesting to have a combinatoric argument which gives 
the form of Theorems 11.21 and 15.11 

Further work will study the fine structure of the asymptotic expansion of the partition function Zn 
when multiple time evolutions are involved. Additionally other partition functions over different families of 
Random Matrices encode similar combinatoric data. For example, in the case when 



where the integral is taken over the space of N x N symmetric matrices, we find that the terms in the 
asymptotic expansion of log (Zn) give generating functions for the number of unoriented maps partitioned 
by the Euler characteristic associated with the embedding of the map. 

What connections these problems will have to the higher Catalan numbers or to other interesting (and 
classical) combinatoric problems remains to be seen. 
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